Charged hadron transverse momentum distributions in Au+Au collisions 

at v /s^ = 200 GeV 



m 
o 
o 

(N 
X> 

PL, 



> 
in 

o 

o 
m 
O 

=5 



B.B.Back 1 , M.D.Baker 2 , D.S.Barton 2 , R.R.Betts 6 , M.Ballintijn 4 , A.A.Bickley 7 , R.Bindel 7 , A.Budzanowski 3 , 
W.Busza 4 , A.Carroll 2 , M.P.Decowski 4 , E.Garcia 6 , N.George 1 ' 2 , K.Gulbrandsen 4 , S.Gushue 2 , C.Halliwell 6 , 
J.Bamblcn 8 , G.A.Beintzelman 2 , C.Henderson 4 , D.J.Hofman 6 , R.S.Hollis 6 , R.Holynski 3 , B.Holzman 2 , A.Iordanova 6 , 
E.Johnson 8 , J.L.Kane 4 , J.Katzy 4 ' 6 , N.Khan 8 , W.Kucewicz 6 , P.Kulinich 4 , C.M.Kuo 5 , W.T.Lin 5 , J.W.Lee 4 , S.Manly 8 , 
D.McLeod 6 , A.C.Mignerey 7 , R.Nouicer 6 , A.Olszewski 3 , R.Pak 2 , I.C.Park 8 , H.Pernegger 4 , C.Reed 4 , L.P.Remsberg 2 , 

M.Reuter 6 , C.Roland 4 , G.Roland 4 , L.Rosenberg 4 , J.Sagerer 6 , P.Sarin 4 , P.Sawicki 3 , W.Skulski 8 , S.G.Steadman 4 , 
P.Steinberg 2 , G.S.F.Stephans 4 , A.Sukhanov 2 , J.-L.Tang 5 , R.Teng 8 , A.Trzupek 3 , C.Vale 4 , G.J.van Nieuwenhuizen 4 , 
R.Verdier 4 , G.I.Veres 4 , BWadsworth 4 , F.L.H.Wolfs 8 , BWosiek 3 , KWozniak 3 , A.H.Wuosmaa 1 , BWyslouch 4 

1 Argonne National Laboratory, Argonne, IL 60439-4843, USA 
2 Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
3 Institute of Nuclear Physics, Krakow, Poland 
4 Massachusetts Institute of Technology, Cambridge, MA 02139-4307, USA 
5 National Central University, Chung-Li, Taiwan 
6 University of Illinois at Chicago, Chicago, IL 60607-7059, USA 
7 University of Maryland, College Park, MD 20742, USA 
8 University of Rochester, Rochester, NY 14627, USA 
(Dated: February 6, 2008) 

We present transverse momentum distributions of charged hadrons produced in Au+Au collisions 
a * ^/ s nn — 200 GeV. The spectra were measured for transverse momenta pr from 0.25 to 4.5 GeV/c 
in a rapidity range of 0.2 < y-n < 1.4. The evolution of the spectra is studied as a function of collision 
centrality, from 65 to 344 participating nucleons. The results are compared to data from proton- 
antiproton collisions and Au+Au collisions at lower RHIC energies. We find a significant change 
of the spectral shape between proton-antiproton and peripheral Au+Au collisions. Comparing 
peripheral to central Au+Au collisions, we find that the yields at high pr exhibit approximate 
scaling with the number of participating nucleons, rather than scaling with the number of binary 
collisions. 
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The yield of charged hadrons produced in collisions of 
gold nuclei at an energy of y / s~vjv = 200 GeV has been 
measured with the PHOBOS detector at the Relativis- 
tic Heavy-Ion Collider (RHIC) at Brookhaven National 
Laboratory. The data are presented as a function of 
collision centrality and transverse momentum px- The 
goal of our measurements is to test our understanding 
of Quantum Chromodynamics (QCD), the fundamental 
theory of the strong interaction. Using nuclear collisions, 
QCD can be studied in a regime of high temperature 
and energy density. Calculations suggest [1| that under 
these conditions a new state of matter, the quark-gluon 
plasma, will be formed. 

In the theoretical analysis of particle production in 
hadronic and nuclear collisions, a distinction is often 
made between the relative contributions from "hard" 
parton-parton scattering processes and "soft" processes. 
The former can be calculated using perturbative QCD, 
whereas the latter are typically treated by phcnomeno- 
logical models that describe the non-perturbative sector 
of QCD 0. The contribution from hard processes 
is expected to grow with increasing collision energy 



and to dominate particle production at high transverse 
momenta. Collisions of heavy nuclei offer ideal conditions 
to test our understanding of this picture, as "hard" 
processes are expected to scale with the number of binary 
nucleon-nucleon collisions N co u, whereas "soft" particle 
production is expected to exhibit scaling with the number 
of participating nucleons N part . In Glauber- model cal- 
culations, N co u scales approximately as (N part ) 4 ^ 3 . For 
central collisions of Au nuclei, this leads to an increase 
of a factor of six in the ratio of N co ii/N part , relative to 
proton-proton collisions. Collision centrality is therefore 
a parameter by which the relative contributions of hard 
and soft processes to particle production can be varied. 
It was found that a description based on a combination of 
hard and soft particle production is compatible with the 
observed centrality dependence of the charged particle 
multiplicity at mid-rapidity in Au+Au collisions 0, Q , 
although the result can also be described by models 
which incorporate initial state parton saturation A 
simultaneous study of the px and centrality dependence 
of particle yields may allow us to better discriminate 
between these different physical pictures. 
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For Au+Au collisions at RHIC energies, it has been 
predicted that the yield and momentum distribution 
of particles produced by hard scattering processes may 
be modified by "jet quenching", i.e. the energy loss 
of high momentum partons in the dense medium 0. 
This phenomenon has been proposed as a diagnostic 
tool for characterizing the parton density in the initial 
stage of high-energy nuclear collisions. First results for 
central Au+Au collisions at y/s NN = 130 GeV have 
shown that particle spectra at large transverse momenta 
indeed change in comparison to p+p collisions Q. This 
result has been used to estimate the energy loss in the 
dense medium However, calculations have not yet 
shown a consistent interpretation of the data on spectral 
shapes and particle multiplicities. The latter seem to 
be overpredicted when the effects of jet quenching 
are included. Here we probe the scaling of particle 
production over a wide range of N co u and N part in 
different regimes of px that are traditionally believed to 
be dominated by soft or hard processes. 

The data were collected using the PHOBOS two-arm 
magnetic spectrometer IS IH The spectrometer 
arms are each equipped with 16 layers of Silicon sensors, 
providing charged particle tracking both outside and 
inside the 2 T field of the PHOBOS magnet. Particles 
within the geometrical acceptance region used in this 
analysis traverse at least 14 of the layers. A two layer 
silicon vertex detector covering \r]\ < 0.92 and 25% of the 
azimuthal angle provided additional information on the 
position of the primary collision vertex. In total, 135168 
detector elements were read out, of which less than 2% 
were non-functional. 

The primary event trigger was provided by two sets of 
16 scintillator paddle counters covering pseudorapidities 
3 < |?7| < 4.5. Additional information for event selection 
was obtained from two zero-degree calorimeters which 
detect spectator neutrons. Details of the event selection 
and centrality determination can be found in |J, Il3| . 
Monte Carlo (MC) simulations of the apparatus were 
based on the HIJING event generator [lj| and the 
GEANT 3.21 simulation package, folding in the signal 
response for scintillator counters and silicon sensors. 

For this analysis, the events were divided into six 
centrality classes, based on the observed signal in the 
paddle counters. Given the monotonic relationship be- 
tween the multiplicity of produced particles in the paddle 
acceptance and the number of participating nucleons 
(Npart), the results of a Glauber-model calculation were 
used 0, 0| to estimate the average number of partici- 
pating nucleons and the number of binary collisions for 
each centrality class. The resulting values are shown in 
Table I. 

As the geometrical layout of the PHOBOS detector 



leads to an asymmetry in the acceptance and detection 
efficiency for positively and negatively charged particles 
for a given magnet polarity, data were taken using 
both polarities. The reproducibility of the absolute 
field strength was found to be better than 1%, based 
on Hall probe measurements for each polarity and the 
comparison of mass distributions for identified particles 
at the two polarities. 

To optimize the precision of the vertex and track 
finding, only events with a reconstructed primary vertex 
position between -10 cm < z vtx < 10 cm along the 
beam axis were selected. By requiring a consistent vertex 
position from the spectrometer and vertex subdetectors, 
in combination with the known position of the beam 
orbit, a vertex resolution better than 0.3 mm (RMS) 
in the longitudinal and vertical directions, z and y, and 
better than 0.5 mm in the horizontal direction, x, was 
achieved. 

Details of the reconstruction algorithm for particle 
tracks can be found in [l5T | . Track candidates were found 
combining a road-following algorithm in the region out- 
side the magnetic field and a transformation algorithm 
for pairs of hits inside the field. The distance of closest 
approach of each reconstructed track with respect to 
the primary vertex (d v t x ) was used to reject background 
particles from decays and secondary interactions. The 
final track selection was based on the % 2 fit probability 
of a full track fit, taking into account multiple scattering 
and energy loss. For further analysis, particles with a 
rapidity of 0.2 < y v < 1.4 were used, assuming the pion 
mass for calculating y^. 

To obtain the invariant yield for charged hadrons, the 
observed transverse momentum distributions were cor- 
rected for the geometrical acceptance of the detector, the 
efficiency of the tracking algorithm and the distortion due 
to binning and momentum resolution. For the acceptance 
and efficiency corrections, the correction factors in px 
and y vary as a function of vertex position and are 
different for the two charge signs in a given field polarity. 
The correction factors were obtained in vertex bins of 
±2.5 cm size by embedding the same MC particle tracks 
into "empty" events and data events. The events were 
reconstructed with the standard reconstruction chain 
with and without the additional embedded tracks. From 
this comparison, a correction factor as a function of pt 
and multiplicity was extracted. Over the centrality range 
covered in this analysis, the reconstruction efficiency 
decreased by approximately 10%. 

MC-simulations based on HIJING and GEANT were 
used to obtain correction factors for the contribution 
of secondary particles to the charged particle spectra. 
The acceptance for secondary and feeddown particles is 
limited to those produced within 10 cm radial distance 
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from the primary collision vertex, as accepted tracks were 
required to have at least one hit in the first two layers 
of the spectrometer. The background contamination is 
further reduced by requiring the particle tracks to have 
d v tx < 3.5 mm. We found that the contamination varies 
from less than 4% at px ~ 0.3 GeV/c to less than 2% 
at pr > 2.5 GeV/c. The corresponding correction was 
applied independent of collision centrality. 

The contamination by wrongly reconstructed 
("ghost") tracks was studied using two methods. 
The reconstruction of HIJING events showed a rate of 
ghost tracks after all quality cuts of less than 4% at 
p T » 0.5 GeV/c and less than 2% at p T » 2.5 GeV/c. 
At higher pr, the available MC statistics did not allow 
a precise determination of the ratio of ghost tracks to 
found tracks. However, MC simulations showed that 
most ghost tracks arise from a wrong match of track 
segments inside and outside the magnetic field. The rate 
of accidental matches was tested by swapping straight 
and curved track segments between the two spectrometer 
arms. We found a contribution of accidental matches of 
less than 1% in the pr > 3 GeV/c region. 

Finally, the spectra were corrected for the effects of the 
finite bin-width and the momentum resolution using an 
iterative procedure. The resolution was determined using 
track embedding and the full reconstruction and fitting 
procedure. For embedded MC tracks we found a mo- 
mentum resolution of Ap/p = 1 — 9% for the momentum 
range of 0.3 to 10 GeV/c. The corresponding correction 
was performed for each centrality bin separately, taking 
into account the non-Gaussian tails in the p^-resolution 
function. At pr ~ 4 GeV/c, the correction varies from a 
factor of 0.92 for the most peripheral events to 0.89 for 
the most central events. 
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Centrality 


~N events 


{Npart} 


(N coU ) 


45-50% 


593942 


65 ±4 


107 


35-45% 


1196402 


93 ±5 


175 


25-35% 


1192357 


138 ±6 


300 


15-25% 


1193373 


200 ±8 


500 


6-15% 


1057007 


276 ±9 


780 


0-6% 


604555 


344 ±11 


1050 



TABLE I: Details of the centrality bins used. The estimated 
uncertainty in (N co ll) ranges from 15% to 10% going from the 
most peripheral to the most central bin. 



In Fig.^we show the invariant yield of charged hadrons 
as a function of transverse momentum, obtained by 
averaging the yields of positive and negative hadrons. 
Data are shown for 6 centrality bins, ranging from 
(Npart) = 65 to {N par t) = 344. The integrated 
yields, when scaled by (N part ), increase by 15% over the 
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FIG. 1: Invariant yields for charged hadrons as a function of 
Pt for 6 centrality bins. For clarity, consecutive bins are scaled 
by factors of 10. Statistical and systematic uncertainties are 
smaller than the symbol size. 



centrality range, consistent with the centrality evolution 
of the mid-rapidity particle density presented in Q . 

In Fig. |21 we show the ratios of data for the most 
peripheral bin ((N par t) = 65 ± 4) and the most central 
bin ((N pa rt) = 344 ± 11) to a fit of the measured pr- 
distribution in proton-antiproton collisions at the same 
energy 16] . In both cases, the Au+Au data were divided 
by the respective value of (Ap art /2). The brackets 
indicate the systematic uncertainty in the Au+Au data 
(90% C.L.). The largest contribution to the systematic 
uncertainty are the overall tracking efficiency, indepen- 
dent of pt, and the px dependent momentum resolution 
and binning correction. The small contamination by 
secondary particles and feeddown particles, due to the 
proximity of the tracking detectors to the collision vertex, 
also ensures that the uncertainty in the corresponding 
correction factors is small. Similarly, the high granularity 
and resolution of the tracking planes leads to small 
uncertainty in the rate of 'ghost' tracks. 

As has been shown previously 0, the yield per 
participant pair in Au+Au collisions at these centralities 
is significantly larger than in proton-antiproton collisions 
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FIG. 2: Ratio of the yield of charged hadrons as a function of 
pr for the most peripheral bin ({N pa rt} = 65 ±4, upper plot) 
and the most central bin ({N par t) = 344 ± 11, lower plot) to 
a fit of proton-antiproton data (see text) scaled by (N par t/2). 
The dashed (solid) line shows the expection of N co u (N par t) 
scaling relative to p + p collisions. The brackets show the 
systematic uncertainty of the Au+Au data. 



at the same energy. We also observe that in peripheral 
Au+Au collisions with (N part ) = 65, corresponding to an 
impact parameter b ~ 10 fm, the spectral shape is already 
strongly modified from that in p+p collisions. It is worth 
noting that the ratio ^ c °' t 'J 2 ) increases by a factor of 
almost three from p + p to the most peripheral Au+Au 
collisions studied here. For the highest pr, the yield for 
central events is significantly smaller than expectations 
based on A^z-scaling. 

The detailed evolution of the spectra from peripheral 
to central events is shown in Fig.|21 for the six centrality 
bins. Each spectrum has been been divided by (N part /2) 
and a fit to the data for the most peripheral bin. N part - 
scaling would correspond to the ratio being constant 
at one (solid line), whereas A^ co ; r scaling, shown by the 
dashed line, would lead to an increase by a factor of two 
from peripheral to central events. Again, the brackets 
show the systematic uncertainty (90% C.L.). It is 
remarkable that the change in spectral shape over this 
range of centralities is small compared to that between 



FIG. 3: Charged hadron yield in Au+Au in six centrality 
bins, divided by a fit to the most peripheral bin and by 
{N P art/2). The dashed (solid) line shows the expectation 
for N co ii (N P art) scaling relative to peripheral collisions. The 
brackets show the systematic uncertainty in the Au+Au data. 



p + p collisions and peripheral events shown in the top 
panel of Fig. |2 In particular at pr > 3 GeV/c, the 
data scale with the number of participants to very good 
approximation. 

The scaling behaviour in various regions of px is 
further illustrated in Fig. Here we show the yield 
per participant pair in six bins of px between 0.45 and 
4.25 GeV/c, normalized to the yield per participant 
pair in the most peripheral bin. The expectation for 
-^Vcozz-scaling relative to the most peripheral bin is again 
shown as a dashed line. The observed evolution can be 
contrasted with the expectation that particle production 
should be characterized by a change from (N part ) scaling 
at low pt to (N co u) scaling at high pt, with (N part ) 
increasing from 65 to 344 and (N co u) from 107 to 1050 
between the most peripheral collisions and the most 
central collisions studied here. No corresponding increase 
in particle production per participant at pr — 3 GeV/c 
and above is observed. Rather, the yields in this region 
scale approximately with the number of participating 
nucleons. 

The failure of binary collision scaling at intermediate 
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FIG. 4: Charged hadron yields per participant pair in 6 
different transverse momentum bins, plotted as a function of 
Np ar t. The data are normalized to the yield in the most 
peripheral centrality bin. The dashed (solid) line shows 
the expectation for N co u (N par t) scaling from peripheral 
to central collisions. The brackets indicate the systematic 
uncertainty for the centrality evolution of this ratio (90% 
C.L.). 



and large px was previously observed in Au+ Au collisions 



a t \/ s nn = 130 GeV Q. This effect has been interpreted 
as the result of the energy loss of jets in the dense 
medium, leading to a suppression of leading hadrons 
relative to the expected scaling. The observed particle 
spectra are believed to depend on a complex interplay 
of initial and final state effects in the production and 
fragmentation of high px partons, including initial state 
multiple scattering, nuclear shadowing and parton energy 
loss. Our observations show that, surprisingly, the com- 
bination of all effects leads to an apparently very simple 
centrality scaling at large pr, providing a challenge to 
theoretical descriptions. That this is not accidental is 
suggested by the agreement with earlier measurements 

at = 130 GeV - Tnc lower ener sy data EH 

when compared over the same centrality range, show very 
similar centrality scaling, even though the invariant yield 
at high pt increases much more rapidly with increasing 
beam energy than the overall particle yield. It has 
recently been argued that the observed scaling could be 
naturally explained in a model assuming the dominance 
of surface emission of high-p^ hadrons |19|. However, the 
approximate participant scaling has also been explained 
in the context of initial state saturation models pcj . 
Experimentally, upcoming studies of d+Au collisions at 
RHIC will provide further insight into the modification 
of particle spectra in a nuclear environment. 
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